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Abstract 

Coincident with shifting monsoon weather patterns over India, the phytoplankter Noctiluca miliaris has 
recently been observed to be dominating phytoplankton blooms in the northeastern Arabian Sea during the 
winter monsoons. Identifying the exact environmental and/or ecological conditions that favor this species has 
been hampered by the lack of concurrent environmental and biological observations on time and space scales 
relevant to ecologic and physiologic processes. We present a bio-optical proxy for N. miliaris measured on highly 
resolved depth scales coincident with hydrographic observations with the goal to identify conducive hydrographic 
conditions for the bloom. The proxy is derived from multichannel excitation chlorophyll a fluorescence and is 
validated with microscopy, pigment composition, and spectral absorption. Phytoplankton populations dominated 
by either diatoms or other dinoflagellates were additionally discerned. N. miliaris populations in full bloom were 
identified offshore in low-nutrient and low-N : P ratio surface waters within a narrow temperature and salinity 
range. These populations transitioned to high-biomass diatom-dominated coastal upwelling populations. A week 
later, the N. miliaris blooms were observed in declining phase, transitioning to very-low-biomass populations of 
non-TV. miliaris dinoflagellates. There were no clear hydrographic conditions uniquely associated with the N. 
miliaris populations, although N. miliaris was not found in the upwelling or extremely oligotrophic waters. 
Taxonomic transitions were not discernible in the spatial structure of the bloom as identified by the ocean color 
Chi imagery, indicating that in situ observations may be necessary to resolve community structure, particularly 
for populations below the surface. 


The Arabian Sea is one of the most productive regions 
in the world because of periods of strong monsoonally 
driven coastal upwelling and convective mixing that bring 
nutrient-rich waters from below the permanent thermocline 
to the surface waters (Wiggert et al. 2005). Monsoons act as 
one of the predominant climatic systems in the Arabian Sea 
and are generated in response to the adjacency of the large 
Eurasian landmass to the Arabian Sea producing intense 
land-sea temperature gradients. The main physical forcings 
that control productivity in the Arabian Sea are winter 
convective mixing during the northeast monsoon (NEM) in 
the winter, December through March, and coastal upwell- 
ing during the southwest monsoon (SWM) in summer, June 
through September (Parab et al. 2006). Winds associated 
with SWM along the northwestern coast of Africa and the 
Arabian Peninsula drive offshore Ekman transport and 
coastal upwelling in the summer and produce large coastal 
phytoplankton blooms (Garrison et al. 1998; Latasa and 
Bidigare 1998; Shalapyonok et al. 2001). The NEM 
produces cooling of the northern Arabian Sea in the 
winter, causing the surface water to sink, generating 
convective mixing, and injecting nutrients to the surface 
(Madhupratap 1999). In effect, these monsoons create two 
distinct blooms, a strong bloom in the summer off the coast 
of Somalia and a weaker bloom in the winter off the 
northwest coast of India and southern Pakistan (Brock 
et al. 1992; Banse and English 2000). 
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The Arabian Sea is a relatively biodiverse ocean basin 
with enhanced biological activity providing large amounts 
of sinking particulate organic matter that maintains an 
extraordinarily extended deep oxygen minimum zone. The 
semi-enclosed basin with convergent and downwelling 
geostrophic circulation characteristic of subtropical gyres 
leads to low oxygen regeneration in the interior deep water 
as well as upwelling driving deep low-oxygen waters to the 
surface around the perimeter (Wajih et al. 2006). 

Recent and significant reductions in snow coverage in 
the Himalayas as well as extremely high summertime 
landmass temperatures are implicated in driving shifts in 
the intensity and timing of sea-level pressure gradients from 
land to sea, intensifying summer monsoons, increasing 
productivity, and further depleting oxygen in the deep 
water of the Arabian Sea (Goes et al. 2005). These oxygen- 
depleted waters are subsequently upwelled to the surface in 
the northeastern Arabian Sea during the winter monsoons, 
and contemporaneous changes in phytoplankton composi- 
tion have been observed (Parab et al. 2006; Gomes 
et al. 2008). The typical pattern in phytoplankton bloom 
succession was dominance by diatoms, followed by cyano- 
bacteria, and then dinoflagellates with the lowest recorded 
concentrations (Sawant and Madhupratap 1996). However, 
since the 1990s the phytoplankton community has become 
dominated by an emergent species, Noctiluca miliaris (also 
known as the green Noctiluca scintillans; Parab et al. 2006; 
Gomes et al. 2008; S. Prakash unpubl.). Although this 
species has been observed extensively in the waters of 
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Southeast Asia, its predominance in the Arabian Sea is a 
recent occurrence (Harrison et al. 2011). 

Combining remotely sensed seasonal bloom distribution 
with shipboard microscopic analysis, Gomes et al. (2008) 
speculate that climatic changes have led to the shift in 
phytoplankton composition from diatoms to dominance 
by dinoflagellates, specifically N. miliaris, in these low- 
oxygenated waters during the winter NEM (February- 
March) and the spring intermonsoon (SIM; Madhu et al. 
2012). Further analysis of the satellite data (Gomes et al. 
2008) reveals that the once-dominant cyanobacteria Tri- 
chodesmium sp. in the eastern Arabian Sea (Westberry and 
Siegel 2006) has been replaced by N. miliaris during the 
SIM and that the N. miliaris blooms are widespread not 
only offshore, but also in the polluted coastal waters off of 
India (Chaghtai and Saifullah 2006). It is thought that 
recent increases in coastal phosphate concentrations give 
N. miliaris a competitive advantage compared to other 
phytoplankton species. Observations of ammonia secretion 
suggest that N. miliaris has sufficient nitrogen for growth 
(Furuya et al. 2006), perhaps because of the combination of 
both heterotrophy and symbiosis. 

N. miliaris — The phytoplankter N. miliaris (Fig. 1) is an 
unpigmented phagotrophic and omnivorous dinoflagellate 
that is large enough to be seen by eye (>500 pm diameter). 
It consumes an array of organisms including other 
phytoplankton, small zooplankton, fish eggs, and detritus. 
N. miliaris consumes its food using a tentacle that connects 
to the cytosome and can obtain food either by direct 
contact with food particles that land on the tentacle or by 
protruding mucus from the tentacle that acts as a net to 
catch food (Kiorboe and Titelman 1998; Chaghtai and 
Saifullah 2006). Although not an active vertical migrator, 
N. miliaris is highly buoyant because of its ability to secrete 
ammonia into the cytoplasm and maintain positive 
buoyancy (Furuya et al. 2006; Lirdwitayaprasit et al. 
2012). Despite being relatively immobile, N. miliaris 
intercepts its food through passive migration (Kahn and 
Swift 1978). The high buoyancy often causes N. miliaris 
blooms to accumulate at the surface, forming noticeably 
discolored waters (Kiorboe and Titelman 1998). 

N. miliaris is able to tolerate moderate fluctuations in 
both salinity and temperature, enabling it to prevail in 
variable hydrographic conditions (Lirdwitayaprasit et al. 
2006). It is suggested that N. miliaris tolerates low oxygen 
conditions because of its symbiont, Pedinomonas noctilu- 
cae , photosynthesizing more efficiently than other phyto- 
plankton in low oxygen conditions (Hansen et al. 2004; 
Gomes et al. 2008). This characteristic facilitates N. miliaris 
dominance during the NEM when oxygen-deficient water is 
upwelled (H. Gomes and J. Goes pers. comm.). 

Classified as a unicellular microscopic flagellate, P. 
noctilucae belongs to the class Prasinophyceae, although it 
gains its green color from the pigments common to the 
Chlorophyceae. Both N. miliaris and P. noctilucae occur 
preferably at low temperatures in tropical waters (Chaghtai 
and Saifullah 2006). The symbiotic relationship exists 
because P. noctilucae ostensibly provides N. miliaris with 
photosynthetic products, which N. miliaris would not have 



Fig. 1. (A) Microscopic image of Noctiluca miliaris collected 

from the Arabian Sea in 2009. (B) Water sample collected from 
the surface waters in the Arabian Sea in 2009 showing the large 
buoyant green N. miliaris cells. Photo credit J. Goes. 

otherwise. In exchange, N. miliaris gives the symbiont 
protection from grazing and provides the low-pH environ- 
ment (4.5-5) critical to P. noetilucae’s survival (Furuya 
et al. 2006). The green pigmentation of P. noctilucae produces 
the distinctive green coloration of N. miliaris so different 
from most other dinoflagellates, specifically the absence of 
peridinin and the presence of chlorophyll b (Chi h). 

Taxonomic identification via pigment-based phytoplank- 
ton functional types — Although chlorophyll a (Chi a) is 
the primary fluorescing pigment, the energy it fluoresces 
originates from light absorption by the suite of accessory 
photosynthetic pigments (Huot and Babin 2010). Fluores- 
cence intensity is proportional to light absorption, via the 
fluorescence quantum yield, and hence to pigment concen- 
tration. By using the excitation and emission characteristics 
associated with phytoplankton pigment fluorescence, the 
composition of phytoplankton communities can be deter- 
mined (Yentsch and Phinney 1985; Lutz et al. 2001; 
Chekalyuk and Hafez 2008, 2013). Because of the distinct 
chlorophyte-like pigmentation, we suspected that N. 
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Fig. 2. Example of calibration slopes (in vivo Chi a 
fluorescence per extracted Chi a concentration), obtained with a 
WET Labs 3X1M and derived as in Fig. 3, for diverse 
phytoplankton species at (A) 440 nm and (B) 532 nm excitations. 
Error bars indicate standard deviation for cultures grown under 
high and low light and collected at different growth phases (data 
modified from Proctor and Roesler 2010). (C) Fluorescence 
pigment ratio (F532 : T 440 ) for each species after transfer function 
applied (see text). The bars for each species is shaded according to 
pigment lineage (black, gray, and white for green, red, and cyano 
lineages, respectively), following Falkowski et al. (2004) and rank 
ordered within lineage according to the fluorescence ratio. 


miliaris would be discernible from diatoms and cyanobac- 
teria based upon not only its pigmentation (Letelier et al. 
1993; Roy et al. 2011), but its associated absorption spectra 
(Hoepffner and Sathyendranath 1993; Sosik and Mitchell 
1995; Johnsen et al. 201 1), and hence its Chi a fluorescence 
excitation (Culver and Perry 1999; Lutz et al. 2001; 
Richardson et al. 2010) in the blue to green range. 

Proctor and Roesler (2010) calibrated the response of 
a multi-excitation (three-channel) Chi fluorometer to 13 
different species of phytoplankton under a range of growth 
phases and irradiance levels. They found that the major 
taxonomic groups could be distinguished from each other 
based upon differential pigment composition. In this paper 
we hypothesize that the pigmentation differences between 
N. miliaris and the diatoms and cyanobacteria previously 
found in the northeastern Arabian Sea are sufficient to 
develop a robust bio-optical proxy based upon multichan- 
nel excitation of Chi fluorescence. Using this proxy, we 
identify the extent of N. miliaris populations in the 
northeastern Arabian Sea in situ during the winter 
monsoon and relate the population distribution to both 
the hydrographic features and the spatial structure of the 
satellite-derived Chi features. We further hypothesize that 
the distributions of N. miliaris will be associated with the 
bloom structures identified by high concentrations of 
satellite-derived Chi in the northeastern Arabian Sea 
during the winter monsoon and that these populations will 
occupy water masses with characteristics distinct from 
other phytoplankton communities. In other words, when 
N. miliaris is observed, it will be in high concentrations and 
will exclude other phytoplankton groups within defined 
water masses because of the competitive advantage that it 
has by both feeding on phytoplankton and receiving 
photosynthate from its symbionts. 

Methods 

Sensor calibration — The Western Environmental Tech- 
nologies (WET) Labs 3X1M is a three-channel (440, 470, 
and 532 nm)-excitation, single-channel (695 nm)-emission 
detection fluorometer in their Environmental Character- 
ization Optics Triplet sensor series. The excitation is 
supplied by low-power light-emitting diodes (LEDs) with 
full-width, half-maximum (FWHM) ranges of 24 nm for 
the 440 nm and 470 nm LEDs and 32 nm for the 532 nm 
LED. The detector optical filter has a FWHM of 70 nm 
and is made of RG645 Schott glass. As the sensor is a 
special-order system with user-defined configurations, it 
was not supplied with calibration factors. Proctor and 
Roesler (2010) quantified the response of the first sensor to 
a variety of phytoplankton species grown under a range of 
controlled growth conditions. Calibration response curves 
(i.e., calibration slope parameters, fluorescence : Chi, units 
digital counts [DC] : [mg m~ 3 ]) varied 5- to 10-fold in 
response to pigmentation differences between species at the 
different wavelengths (Fig. 2A,B). The median calibration 
value for the 470 nm excitation (the typical excitation 
wavelength for most commercially available Chi fluorom- 
eters) coincided with those obtained for Thalassiosira 
pseudonana, a ubiquitous centric diatom in the red pigment 
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Fig. 3. Standard calibration curves for each excitation 
channel of the WET Labs 3X1M sensor using a dilution series 
of Thalassiosira pseudonana. The y error bars represent the 
standard deviation of a 60 s burst sample of raw fluorescence for 
each dilution, and the x error bars represent the standard 
deviation of triplicate extracted Chi a analyses for each dilution. 
Calibration slopes are computed from type II linear regression. 

lineage. Dinoflagellates typically had higher calibration 
slopes and cyanobacteria and chlorophytes had lower 
slopes. In addition, the calibration slopes obtained for T. 
pseudonana were repeatable within 5% among calibrations 
separated by years as long as the growth conditions were 
replicated. This provides a useful means for computing a 
transfer function between sensors, essentially adjusting for 
difference in gain at each excitation wavelength. After 
adjustment by the transfer function, the species-specific 
calibration factors obtained by Proctor and Roesler (2010) 
can be computed for any sensor with the same optical 
configuration. We used this approach to recalibrate both 
the original sensor (serial number SN 001) and a newer 
sensor (SN 004), the latter of which was deployed in the 
Arabian Sea. 

An acclimated culture of T. pseudonana (CCMP 1335) 
was grown in sterile filtered LI medium (a vitamin-, trace 
metal-, and nutrient-enriched seawater; Guillard and 
Hargraves 1993, as modified by Guillard and Morton 
2003) under controlled light conditions (12:12 light: dark 
cycle at 250 /(mol photons m -2 s -1 ) and harvested in 
exponential growth at a concentration of approximately 
20 mg m -3 Chi a. A series of 10 dilutions of the culture was 
made, starting with 100% culture and decreasing by a 
factor of two via dilution with the same filtered seawater 
used for media. Two blanks were also run: Barnstad 
Nanopure® water and the filtered seawater diluent. 
Subsamples from each of the 10 culture dilutions were 
filtered for triplicate extractive Chi a analysis (Yentsch and 
Menzel 1963). 

The absorption spectrum for the culture of T. pseudo- 
nana was determined spectrophotometrically from a 
subsample of the culture filtered onto a glass fiber filter 
(Whatman GF/F). The particulate absorption analysis was 


performed on the Cary 300 Ultraviolet visible spectro- 
photometer following Kishino et al. (1985) for separation 
of phytoplankton and non-algal particle absorption and 
Roesler (1998) for configuration and absorption calcula- 
tions. The absorption spectrum for T. pseudonana was 
then compared to that for previous calibrations of 
instrument SN 001 (Proctor and Roesler 2010) to ensure 
consistent absorption ratios (within 5% coefficient of 
variation) at the fluorometer excitation wavelengths 
between calibrations. 

A type II regression analysis of the median fluorescence 
values for each wavelength at 440 nm, 470 nm, and 532 nm 
and the extracted Chi a concentrations for each dilution 
was performed (Fig. 3). The linear relationship is of the 
form Y = M X X + B, where Y is the observed fluorescence 
(units DC), X is the concentration of the calibration 
component Chi (mg m -3 ), M is the slope of the response 
(DC : [mg m 3 ] ), and B is the y-intercept, representing the 
signal from the media blank (DC). Instrumental dark 
values, F dark (units DC), were measured with the optical 
face of the sensor covered in electrical tape and submerged 
in the media as recommended by the manufacturer. From 
this calibration, in situ Chi a concentration can be 
estimated from the sample fluorescence, F sample , using the 
expression 

Chi a (mg m -3 ) = (^sample — ^dark )/^ ( 1 ) 

where F dark is the dark reading of the sensor at the specific 
wavelength and M cal is the calibration slope determined for 
the T. pseudonana culture for that specific wavelength. 
Each time the sensor is used to quantify Chi, a media blank 
reading is collected using filtered sample (either from a 
culture or from the natural environment). The reading is 
processed as if it were a sample, the media-specific Chi 
value, which is subsequently removed from each sample 
value (Davis and Cullen 2003). Typically this value is not 
significantly different from zero; however, exceptions are 
found, particularly in waters with high-colored dissolved 
organic matter (Proctor and Roesler 2010). 

Transfer functions between sensors — Just as sensors are 
calibrated through time to account for instrument drift, so 
too are multiple sensors intercalibrated with the same 
standards to quantify and account for differences in gains 
(response curves) for each excitation/emission pair. By 
computing the ratio of the gains for the T. pseudonana 
culture, a transfer function between the two 3X1M sensors 
(SN 001 and SN 004) was computed at each wavelength 

(Tj- 2 ): 

n~ 2 = M '/ M i ( 2 ) 

where M) and M? are the T. pseudonana calibration slopes, 
M ca i, for sensors 1 and 2 measured at the same excitation 
wavelength X. The transfer function is then applied to 
measured fluorescence observations of sensor 2 at wave- 
length X, F^ eaS/l , so that the species calibration table of 
Proctor and Roesler (2010) can be used to compute the 
corrected fluorescence for sensor 2 at wavelength X, F c 2 orr; : 
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Fig. 4. Seven-day composite images of (A,B) Chi a and (C,D) sea surface temperature (SST) from the NASA Moderate Resolution 
Imaging Spectroradiometer onboard the Aqua spacecraft (MODISA). We show two temporal ranges that combine to encompass the field 
campaign: (A,C) 07-13 March 2011, encompassing Sta. 1 to 6, and (B,D) 13-19 March 2011, encompassing Sta. 7 to 13. We binned all 
data to 2 km spatial resolution. We considered all valid MODISA Chi a and SST retrievals when generating these composite images. 
Black circles indicate in situ sampling stations. Regions of gray and black indicate land and no satellite retrievals, respectively. We 
acquired these MODISA data from the NASA Ocean Biology Processing Group (OBPG) at Goddard Space Flight Center (http:// 
oceancolor.gsfc.nasa.gov/cgi/browse. pi) and spatially and temporally binned the imagery using standard NASA OBPG software (http:// 
oceancolor.gsfc.nasa.gov/seadas/). 


F 2 =F 2 xT 1 " 2 (3) 

corr/l meas/l X V / 

The transfer function serves to correct for differences 
between sensors and also for a single sensor between 
calibrations if the gains are drifting. As new species are 
included in the calibration table it is necessary to tie them 
to a calibration with T. pseudonana or another reference 
species already in the table. 

Although the fluorescence response ( M for each species) 
varies between species depending upon excitation wave- 
length, the fluorescence ratios are statistically different 
between the pigment lineages (Fig. 2C), and in some cases 
within lineages. The calibration slopes often do not group 
by lineage (Fig. 2A,B), but the ratios clearly do. The 
calibrated fluorescence ratio F 532 '■ F 440 clearly exhibits very 
low values for the green lineage relative to the red lineage 
of the diatoms and dinoflagellates, whereas that of the 
cyanophyte lineage exhibits very high values compared to 
the red lineage. Given that the fluorescence ratios include 
such within-species pigment ratio variations associated 
with photoacclimation and growth phase (shown by error 


bars), we predict that the fluorescence ratios associated 
with N. miliaris will be most similar to the green lineage and 
therefore significantly different from the red and cyano 
lineages. 

Field observations — In situ profiles of hydrographic and 
optical observations as well as water samples were collected 
at 13 stations aboard the R/V Sugar Sampada (SS286) from 
06 to 21 March 2011 in the northeastern Arabian Sea 
(Fig. 4). At each station, fluorescence profiles were 
measured in concert with a Falmouth Scientific Instru- 
ments Micro conductivity, temperature, depth sensor. 
Discrete water samples were collected via a Niskin rosette 
for analyses of nutrients, pigments, absorption, and 
microscopy. 

Nitrate, nitrite, phosphate, and silicate were determined 
by autoanalysis. Pigment composition was determined by 
high-performance liquid chromatography (HPLC) by the 
National Aeronautics and Space Administration’s Field 
Support Group. Spectrophotometric absorption of partic- 
ulates was measured using the quantitative filter technique 
(Mitchell 1990) as modified by Roesler (1998). Phytoplank- 
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Fig. 5. Depth contour of (A) temperature, (B) salinity, and (C) Chi a concentration 
(calibrated fluorescence) measured at 13 stations along the cruise track. Maximal Chi of 
~4.5 mg m -3 was observed at Sta. 7. 

ton and non-algal fractions were determined by extraction 
(Kishino et al. 1985). Phytoplankton species were identified 
and quantified microscopically (Parab et al. 2006). 

Results 

The cruise began offshore in stratified waters with 
temperature ranging from 23.5°C to 26.8°C and salinity 
ranging from 35.2 to 36.3 (Fig. 5). This region appeared to 
be a frontal zone with stations drifting between warmer 
surface waters at Sta. 3 and cooler surface water at Sta. 1 
and strongest thermal gradients at Stas. 2, 4, and 6. The 
highest salinities were observed at depth of Stas. 1 and 3 
with lower salinities observed at depth of Stas. 2 and 4-6. 

Onshore at Sta. 7, the water was well mixed and exhibited 
the lowest temperature and salinity values because of 
strong upwelling. The return to the offshore stations 
showed overall surface warming with the highest temper- 


ature at Stas. 1 1 and 12 (Fig. 4C,D). Higher surface salinity 
(like Stas. 2 and 4-6) was found at Stas. 8, 12, and 13, 
whereas Stas. 9-1 1 exhibited higher salinity throughout the 
water column with some evidence of slightly freshened 
conditions at 40-50 m. Hydrographic conditions at Sta. 8 
were most similar to Sta. 6, whereas Stas. 9-1 1 were most 
similar to Stas. 1, 3, and 4; Stas. 12 and 13 were 
hydrographically more similar to Stas. 2, 5, and 6. 

An offshore surface phytoplankton bloom was observed 
at Stas. 1-5 (Fig. 5C), with highest concentrations at Sta. 2 
and lowest at Sta. 3. A subsurface population was observed 
at about 30-50 m at Stas. 3-5 and appeared hydrographi- 
cally connected to the large coastal bloom at Stas. 6 and 7. 
There was a slight subsurface increase in Chi at Sta. 1 as 
well, but not as defined; the deep population was absent at 
Sta. 2. Upon return to the site of the initial surface bloom 
(Stas. 8-10) there was a distinct absence of the previous 
bloom; only a small deep population at approximately 45 m 
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Fig. 6. Depth distribution of discrete (A) nitrate, (B) silicate, and (C) phosphate 
concentrations (/(mol L >) measured at 13 stations along the cruise track. (A) Black line 
indicates the approximate top of the nitricline as determined by the maximal concentration 
gradient. (B,C) Contours indicate three molar ratio values of nitrate to either silicate (0.5, 1.0, and 
2.0) or phosphate (1, 10, 16) concentrations; dashed, solid, and bold indicate lowest to highest 
molar ratio values. 


at Sta. 8 was observed. The transect south evidenced two 
discrete populations, one at 60 m at Sta. 12 and a lower 
concentration at 35 m at Sta. 13. These populations were 
not observed at neighboring stations but were observed 
with two different Chi fluorometers, indicating that they 
were likely not artifacts. 

Nutrient concentrations were generally low, with the 
depth of the nutricline mostly associated with the depth of 
the 25°C isotherm (Fig. 6). Both silicate and phosphate 
were generally higher in the offshore stations the first week, 
decreasing by week 2. The surface waters were relatively 
depleted in inorganic nitrogen relative to both silicate and 
phosphate; the latter was not observed to approach the 
Redfield molar ratio of 16 within the upper 60 m (Fig. 6C). 

The fluorescence ratios provide insight into the phyto- 
plankton community composition (Fig. 7). Although there 
is some spatial coherence between each ratio and some 
features in the Chi distribution, neither ratio exhibits 


coherence over the whole transect. The ratio of F 410 : /’ 440 
generally increased with depth and ranged in value from 
about 0.8 to 1. The ratio of F 532 \F 440 had a much broader 
range, generally 0.25-0.8, although some higher values 
were observed in the lower biomass waters. The ratios for 
the near-surface Chi patch in Stas. 1 and 2 were distinct 
from those observed in deeper populations and were more 
similar to the values offshore of the coastal upwelling 
stations. The mid-depths of Stas. 3-5 had a population that 
differed significantly from the bordering population above 
and below it. The F 532 ■ F 440 values were consistently much 
larger and most similar to that found at depth at Sta. 2. 
Discrete populations were recognizable in the small patches 
as observed at the bottom of Stas. 8 and 12 and at the 
middle depths of Sta. 13. 

Distribution of phytoplankton groups — Profile fluores- 
cence ratios in 2 m bins were plotted in ratio : ratio space 
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Fig. 7. Depth contour of the Chi fluorescence ratios (A) Fn 10 : F 44 0 and (B) F 332 '■ F 44 o- White 
circles indicate discrete water sampling depths used for proxy validations. 


(F 53 2 : F 44 q vs. /* 47 o : F 440 ) to look for characteristic group- 
ings based on taxonomic dominance (Fig. 8). Near the 
surface (<15 m), there were three distinct clusters of data 
points. The ratio values for those three clusters are 
consistent with those observed in laboratory cultures for 
diatoms and dinoflagellates in the red lineage and for the 



Fig. 8. Relationship between F 470 : F 44{) and F 532 : /Aro for 
each 2 m depth bin along the 13 stations. Symbol color represents 
depth (m). The median ratios computed for discrete depths 
coincident with numerical dominance by each phytoplankton 
group are indicated by black symbols; error bars indicate standard 
deviations (values from Table 2). 


green lineage. From those three clusters of surface water 
observations, the values of both fluorescence ratios increase 
nearly monotonically as depth increases. Because both 
ratios have F 440 in the denominator (essentially scaling the 
accessory-pigment-excited fluorescence by Chi biomass, 
which dominates the 440 nm signature), this pattern 
indicates an increase in accessory pigments relative to Chi 
a as depth increases or light decreases. This suggests that 
the ratios may then become proxies for photoacclimation. 
A cluster analysis of the fluorescence ratios and depth 
reveals that 89% of the variability is explained by variations 
in F 532 :F 440 , 8% by variations in F 470 :F 440 , and the 
remainder by depth variations. The three clusters of points 
were statistically different at p < 0.001, indicating that the 
clusters are distinct populations with respect to their paired 
fluorescence ratios. 

Identification — In order to identify which populations 
were associated with the three clusters identified in 
ratio : ratio space, the ratios were examined for those 
discrete depths at which absorption, pigments, and 
microscopic evaluations were made. When N. miliaris 
represented greater than 30% of the total phytoplankton 
community (by number), the fluorescence ratio F 532 ■ F 440 
converged to a characteristic value of 0.60 ± 0.03 
(Fig. 9A). Even when N. miliaris represented only greater 
than 10% of the community by number, the value was 
0.56 ± 0.06. Similar results were found for diatom- 
and dinoflagellate-dominated communities for each ratio 
(Table 1). The ratios were statistically distinct between the 
dominant three groups and, when plotted on top of the in 
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Fig. 9. (A) Chi fluorescence ratio F 532 : F 440 and (B) HPLC-determined pigment ratio of Chi b : Chi a as a function of percentage 

N. miliaris by cell number. 


situ 2 m binned observations, coincide with the three 
clusters (Fig. 8). 

In order to link the species composition to the 
fluorescence ratios via pigmentation, characteristic pig- 
mentation ratios were examined. Chi b in P. noctilucae 
symbionts was expected to provide the strongest pigment 
indicator of N. miliaris. A comparison between the 
percentage N. miliaris by number and the Chi b : Chi a 
ratio shows that there is a strong linear relationship 
between Chi b and N. miliaris (Fig. 9B). The Chi b to Chi 
a ratio observed for N. miliaris ranged from approximately 
0.05 to 0.2 (mg Chi b: mg Chi a) and appeared as two 
different groupings with similar slope but different 
intercept values (approximately 0.08 and 0.14 [mg Chi 
h : mg Chi a\). These roughly correspond to surface and 
deep stations (i.e., high-light and low-light conditions, 
respectively). However, presence of N. miliaris alone is not 
sufficient to explain Chi b concentrations. In fact, there 
were seven samples with ratios above 0.15 that had no N. 
miliaris, indicating the presence of other Chi b containing 
taxa. Interestingly, this other Chi b containing community 
did not have the same fluorescence ratio signature as N. 
miliaris and thus would not be misidentified optically. 

The discrete samples were separated according to 
microscopic dominance by N. miliaris, diatoms, and 
dinoflagellates and the associated phytoplankton absorp- 

Table 1. Characteristic median and standard deviation 
fluorescence ratios for the three major taxonomic groups 
observed at depths coincident with discrete samples analyzed 
microscopically. Noctiluca miliaris values computed for cases 
when dominance exceeded 30% by cell number (Fig. 9A). Diatom 
and dinoflagellate values computed for cases when dominance 
exceeded 60% of population and sample depths were shallower 
than 30 m. 


Phytoplankton group 

F532 : F440 

7*470 ■ 7*440 

N. miliaris 

0.60±0.03 

1.10±0.02 

Diatom 

0.42±0.13 

1.12±0.05 

Dinoflagellate 

0.94±0.22 

1.08+0.06 


tion spectra distinguished by dominant group (Fig. 10). 
The absorption spectra cluster significantly according to 
microscopic evaluation, with N. miliaris- dominated sam- 
ples exhibiting the highest blue to red absorption peak ratio 
and the steepest blue to green slope (e.g., 440-532 nm), 
followed by the dinoflagellates, then the diatoms. The 
systematic decrease in blue to green absorption slope is 
comparable to a systematic increase in green to blue 
absorption or accessory to Chi a pigment absorption. Thus 
the clustering observed in the fluorescence ratios is 
consistent with the clustering observed in absorption 
spectra. 

Application — Starting with the characteristic F 470 :F 440 
and F 532 'F 440 values determined for the samples with 
nearly pure community composition (Table 1) and extend- 
ing those values to the depth_dependent clusters identified 
in ratio : ratio space (three clusters identified in Fig. 8), 
each 2 m bin of the in situ fluorescence profiles was 
assigned a dominant taxonomic group (Fig. 11). A general 
characterization of the stations sampled in the Arabian Sea 
as a function of depth and station indicated that there was 
an evident species shift from the beginning offshore stations 
to the ending offshore stations. N. miliaris was prevalent at 
the surface and mid-depths of all beginning offshore Stas. 
1-5, except Sta. 2, whereas diatoms primarily characterized 
the deep populations of those stations (except at Stas. 1 and 
2) and dominated the high-biomass coastal upwelling Stas. 
6 and 7. The return to the offshore bloom (Stas. 8-10) 
yielded warm and salty stratified waters, quite different 
conditions of very low Chi concentrations, and a commu- 
nity dominated by dinoflagellates. The very same Chi and 
community characteristics were observed at Sta. 2 the week 
before. A remnant of the N. miliaris underlain by diatoms 
was observed below 40 m at Sta. 8. 

Fluorometrically, cyanobacteria were identified in a few 
locations; however, some species are too small to identify 
microscopically for direct validation. An examination of 
the HPLC-derived pigment composition in those regions 
indicated the presence of the pigment zeaxanthin, poten- 
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Fig. 10. Phytoplankton absorption spectra scaled to Chi a absorption peak at 676 nm. 
Spectra are color coded based upon dominance by taxonomic group determined microscopically: 
diatoms (blue), dinoflagellates (red), and N. miliaris (green). Average spectra for each group are 
shown in black. 


tially indicating the existence of phytoplankton of the blue- 
green pigment class (most likely cyanobacteria) at these 
specific locations. However, biomass was very low in these 
waters, and thus the signal to noise ratio was larger. 
Although zeaxanthin is not unique to cyanobacteria and 
could therefore be indicative of other species, the high 
fluorescence ratios support some cyanobacterial presence. 

Validation — Modeled fluorescence data were compared 
to measured microscopy data to confirm successful 
taxonomic identification at the discrete sample locations 
(Table 2). The fluorescence ratios falsely indicated the 
dominance of N. miliaris on two occasions; in one of those 
samples N. miliaris was a large portion of the community 
numerically, and, because of its large size, likely was 


optically dominant. In no case was a population of N. 
miliaris present and not identified by the fluorescence 
signatures. Fluorescence distinguished dominant taxonom- 
ic group 80% of the time, and of the 20% for which 
fluorescence-based dominance did not match microscopi- 
cally measured dominance, it occurred either in very-low- 
biomass areas emitting limited fluorescence signal, or, in 7 
of the 12 cases, there was a significant contribution by the 
group identified fluorometrically. It is notable that the 
diatoms are distinguishable from dinoflagellates by this 
approach despite somewhat similar absorption spectra, a 
result consistent with those of Proctor and Roesler (2010). 
Only eight of the discrete water samples identified N. 
miliaris microscopically, but fluorescence identified contin- 
uous populations, most of which were missed by the low- 



Fig. 11. Spatial distribution and dominance of N. miliaris and other taxonomic groups as predicted by the paired Chi fluorescence 
ratios. Each symbol represents a 2 m depth bin. 
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Table 2. Validation of modeled vs. measured taxonomic 
identification. Comparison between the modeled phytoplankton 
group estimated from the fluorescence ratio proxy and the 
measured phytoplankton group determined by microscopy. 
Numbers in bold indicate the number of groups identified that 
matched the modeled and measured data. Numbers in parentheses 
indicate number of cases for which the measured community had 
significant contributions by the predicted dominant group. 


Modeled : measured 

N. miliaris 

Diatom 

Dinoflagellate 

N. miliaris 

6 

1(1) 

1 

Diatom 

0 

27 

0 

Dinoflagellate 

0 

10(6) 

16 


resolution water sampling, such as the deep populations at 
Sta. 3 and 5. 

Distribution of N. miliaris — The distribution for N. 
miliaris predicted by the fluorescence ratio proxy suggests 
some continuous hydrographic region along Stas. 1, 3-5, 
which is tightly linked to distribution of diatoms in the 
higher biomass regions but transitions to dinoflagellates as 
biomass decreases. However, these features are not so 
clearly outlined in the satellite-derived Chi image, as Sta. 2 
does not appear to be significantly different from Stas. 1 or 
3 whereas Sta. 5 does appear to be a separate feature. The 
return track exhibits large variations in the Chi features from 
Sta. 8 to Sta. 13, which are lower Chi concentrations but 
appear to be dominated by a single group (dinoflagellates). 

Hydrographically, N. miliaris are not uniquely found in 
waters of distinct temperature or salinity ranges (Fig. 12) 
although the observed ranges are narrower than for the 
other detected groups. With the exception of two data 
points, the temperature range is about 24. 25-25. 5°C and 
the salinity range is 35.75-36.25. But within these same 
temperature and salinity ranges the other three plankton 
groups are observed to be dominant at some locations 
within the cruise transect. 



Salinity 


Fig. 12. Temperature-salinity diagram for the 2 m bin 
averaged values observed along the cruise transect. Symbols are 
coded by dominant phytoplankton groups as in Fig. 1 1 . 



PCI 

Fig. 13. Cluster analysis of discrete observations of temper- 
ature, salinity, nutrients, and log-transformed microscopic cell 
counts grouped into diatoms, non-Noctiluca dinoflagellates, and 
Noctiluca. Eighty-three percent of the variation is explained by the 
first two components. 

A cluster analysis of the discrete observations of 
nutrients, hydrography (temperature and salinity), and 
microscopic cell concentrations of the three phytoplankton 
groups (log transformed) indicates that 83% of the 
variability in the distribution of the phytoplankton groups 
can be explained by two components (Fig. 13). Diatoms 
appear to be more covariant with colder waters that have 
higher nutrient concentrations, whereas N. miliaris and 
other dinoflagellates appear to covary and prefer the higher 
salinity waters that are characterized by low N : P ratios. 

Discussion 

A three-wavelength excitation Chi fluorometer, of a type 
that was previously calibrated to quantify the Chl-specific 
fluorescence response (i.e., the calibration slope) for 13 
species in the laboratory (Proctor and Roesler 2010), was 
deployed in the Arabian Sea with the goal to develop a 
fluorescence-ratio-based optical proxy for the elusive 
dinoflagellate N. miliaris. Isolation of the Arabian Sea N. 
miliaris strain is challenging (S. Parab pers. comm.) and the 
symbiotic relationships observed in the Arabian Sea and 
adjacent regions are likely to vary tremendously depending 
upon environmental conditions (Saito et al. 2006; Sriwoon 
et al. 2008). Considering these constraints we used an in 
situ approach to develop the fluorescence ratio proxy for N. 
miliaris, statistically determining the ratio from sample 
depths at which it numerically dominated the phytoplank- 
ton community. Because of its exceedingly large size (and 
hence absorption cross section) the fluorescence proxy 
values were statistically significant above a numerical 
dominance of 10% of the total population. 

A significant outcome of this analysis was the capability 
of using this approach to statistically distinguish diatoms 
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from dinoflagellates, although because of the overall 
smaller cell diameters and absorption cross sections, the 
proxy values were determined when each group was found 
in a higher percentage of numerical dominance. Although it 
has been previously demonstrated that certain species of 
diatoms and dinoflagellates were distinguishable under 
laboratory culture conditions (Proctor and Roesler 2010), 
extending these results to the field with mixed communities 
and a range of environmental conditions is encouraging. 
The pigment differences between the green and red lineage 
phytoplankton groups is significant, whereas the predom- 
inant pigmentation differences between diatoms and 
dinoflagellates are in their carotenoids (fucoxanthin and 
peridinin), which have somewhat similar absorption spectra 
(Roy et al. 2011). 

The strength of the fluorescence ratio optical proxy 
approach is the capability for direct observation of the 
distribution of target phytoplankton species and/or pigment- 
based phytoplankton functional types (PFTs) measured on 
time and space scales consistent with hydrographic obser- 
vations and on scales relevant to the ecological and 
physiological process of the phytoplankton. As target 
species and communities transition between dominance by 
one group over another, the continuous evolution of the 
proxy values from one to another provides the context for 
single observations. In situ validation improves confidence 
in the proxy signal. 

The capability for discerning the depth distribution of 
phytoplankton groups via an optical proxy yields an 
approach for rapidly assessing phytoplankton populations 
from remote platforms such as moorings, floats, and 
gliders, among others, and can provide the necessary 
validation data sets for the next generation of ocean color 
sensors, such as the hyperspectral radiometer that is 
expected to comprise the National Aeronautics and Space 
Agency’s Pre-Aerosol Clouds and Ocean Ecosystem 
mission (decadal.gsfc.nasa.gov/pace.html). Hyperspectral 
remote sensing promises to improve capabilities for 
determining pigment-based PFTs but requires depth 
distributions of phytoplankton populations to validate 
inversion algorithms (Werdell et al. 2014). Such validation 
data sets have been difficult to match up to satellite 
observations historically (Werdell and Bailey 2005) but are 
progressing with improved bio-optical observations and 
proxies for phytoplankton community structure (Werdell 
et al. 2013). 

The continuous profiles of the fluorescence ratios in the 
Arabian Sea were instrumental in revealing information 
about the distribution and dynamics of the species N. 
miliaris with respect to the phytoplankton community. The 
approach allowed us to separate phytoplankton groups in 
situ as a function of depth. At the same time, we observed 
that within the range for a specific species or phytoplank- 
ton group, both fluorescence ratios increased as depth 
increased, consistent with the photoacclimation strategy of 
preferential increase in accessory pigments relative to Chi a 
as incident irradiance decreased. That the taxonomic 
distinctions remained separated as photoacclimation oc- 
curred indicates the robustness of the approach and the 
ability to use it as a means of assessing acclimation history 


temporally. Other environmental conditions that affect 
pigment composition, such as nutrient limitation and 
starvation (Henriksen et al. 2002), are also likely to yield 
variations in proxy values and by inversion provide additional 
physiological information under certain conditions. 

Contrary to our hypotheses, the populations of N. 
miliaris were not uniquely associated with a particular 
water mass or water characteristics. The temperature and 
salinity ranges in which N. miliaris was found were 
narrower than but encompassed the ranges identified for 
phytoplankton communities dominated by diatoms and 
dinoflagellates, with the exception of the strong upwelling 
conditions at Sta. 7. However, it was determined that 
gradients in phytoplankton composition somewhat corre- 
sponded with biomass gradients — that is, as the ship moved 
out of the original blooms and later returned to its 
remnants, there was a decrease in biomass as well as a 
shift in the dominant phytoplankton community and as 
biomass decreased with depth, there was often a change in 
phytoplankton dominance. 

The R/V Sagar Sampada cruise began by entering into a 
large N. miliaris bloom that was predominantly at the 
surface and mid-depth waters. These waters are typically 
oligotrophic and maximal Chi concentrations were 
<2 mg m -3 . A large bloom of diatoms appeared about 
300 km from shore at a depth of 40-50 m at Stas. 3, 4, 
and 5 and Sta. 8 a week later. These offshore diatom 
populations were associated with the N. miliaris popula- 
tions but separated by depth. This diatom population 
appeared contiguous with the large bloom in the upwelling 
along the Gujarat coast. 

The exception to this pattern was observed at Sta. 2, a 
location that was identified at the time as the ship drifting 
out of the bloom. This station was characterized by very 
different hydrography compared to Stas. 1 and 3, and a 
small surface diatom population was underlain by very low 
Chi concentrations dominated by dinoflagellates. From 
satellite this station was shown to be located at the edge of 
the bloom, and the conditions of Sta. 2 were remarkably 
similar to those observed at Stas. 4, 6, and 8, albeit with 
much lower Chi concentrations and very different phyto- 
plankton composition. A week later the remaining tendrils 
of enhanced satellite-detected Chi were associated with low 
concentrations of dinoflagellates and the vestige of a deep 
diatom population likely extending from the coastal 
upwelling bloom. Because the phytoplankton composition 
transitioned from N. miliaris to other dinoflagellates upon 
returning to the offshore stations, we suggest that the N. 
miliaris bloom had terminated by this time, allowing the 
prevalence of other species to occupy the water column. 
The presence of salps, with N. miliaris cells observed in 
their guts, confirmed that the original N. miliaris bloom 
had evolved to a demise stage. 

Satellite observations confirm a declining bloom feature 
from the first to the second weeks of the field campaign as 
well as significant warming of the surface waters. However, 
there is nothing in the remotely sensed Chi features that 
would indicate the transitions in phytoplankton dominance. 
An analysis of these features using semi-analytical inversion 
of remotely sensed reflectance suggested that N. miliaris 
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populations appear on the fringes of diatom blooms (P. J. 
Werdell unpubl.), features that are easily missed by 
traditional shipboard sampling, but were captured in the 
continuous profiles of multichannel fluorescence. That N. 
miliaris was observed in the transition zones between well- 
developed diatom blooms and the low concentrations of 
dinoflagellate populations in the oligotrophic regions 
suggests a sequencing of the three groups that plays out 
both temporally at a single location and spatially at any 
given time. N. miliaris may be at once both grazing upon the 
declining diatom bloom and expanding its distribution in 
low N : P waters by secreting ammonia into the cytoplasm to 
alleviate nitrogen stress in its symbionts. Both processes are 
supported by the observations here. 

Thus three distinct ecological regions are identified: (1) 
the initial N. miliaris bloom associated with satellite- 
derived Chi features, located in cooler and saltier waters 
to about 40 m; (2) diatom populations associated with the 
colder and fresher coastal upwelling features and the 
offshore extensions of those blooms, which are adjacent 
to the N. miliaris bloom initially and the dinoflagellate 
populations in later stages; and (3) the dinoflagellate- 
dominated community inhabiting the fringes of the 
satellite-derived Chi features and into the increasingly 
stratified, oligotrophic waters characterized by very warm 
and salty waters. The latter populations are found in the 
same salinity ranges as N. miliaris but in warmer 
temperatures. In this region the warmer temperatures are 
likely proxies for age of the bloom in stratified waters 
as temperatures of the surface increase with time. Thus, 
the transitions between phytoplankton communities are 
not necessarily associated with hydrographic features 
or processes such as mixing and stratification, but are 
indicative of the temporal evolution of species succession. 

That said, there are limitations to identifying changes in 
N. miliaris dominance both through the fluorescence proxy 
and through satellite imagery. N. miliaris is no longer 
identifiable once it sheds its symbionts, having lost its 
pigments and therefore its fluorescence proxy. N. miliaris 
identification by satellite is limited to the upper depths of 
the water column, emphasizing the need for in situ 
validation such as fluorescence proxies. Although N. 
miliaris is not an active vertical migrator, its ability to 
secrete ammonia and maintain positive buoyancy prolongs 
its exposure to surface water while it is in its growth phase. 
Satellite imagery provides insights into phytoplankton 
dynamics that cannot be determined by traditional in situ 
sampling methods and indicates specific places and times 
for target species activity (Werdell et al. 2014; P. J. Werdell 
unpubl.). Information obtained from satellite is used in 
conjunction with the fluorescence proxy to understand 
widespread trends of N. miliaris distribution. 

Although the oxygen data are not available, comparing 
the oxygen levels between the beginning and ending 
offshore stations would be significant in considering 
whether N. miliaris inhabits low-dissolved-oxygen waters 
and/or exacerbates low-oxygen zones by creating large 
areas of decomposition and subsequent oxygen utilization 
at depth during its decay. Because N. miliaris’’ supportive 
environment continues to widen because of increasing 


oxygen minimum zones as well as the expanding zone 
affected by agricultural runoff from the Indian subconti- 
nent (D’Silva et al. 2012), having continuous observation of 
N. miliaris will be extremely important in the immediate 
and near future in order to best understand the anthropo- 
genic and natural mechanisms supporting this species. 

Contrary to our hypothesis, the populations of N. 
miliaris did not necessarily correspond with high-biomass 
features detected from satellite. Although the highest 
biomass was typically associated with the diatom popula- 
tion (driven largely but not exclusively by the coastal 
upwelling stations), the next highest biomass regions were 
associated with N. miliaris and the lowest with dinoflagel- 
lates. There was not, however, a clear demarcation in 
phytoplankton group dominance that corresponded exclu- 
sively to the Chi features identified in the remotely sensed 
imagery, suggesting that phytoplankton composition in this 
region did not strictly control total biomass. In other 
words, biomass was not a robust marker of composition as 
has been found for other species such as Karenia brevis in 
the waters around Florida (Stumpf et al. 2003). This 
implies that remote sensing approaches to detect N. miliaris 
cannot rely on biomass indicators but must incorporate 
spectral signatures associated with pigment absorption 
features and perhaps reduced backscattering efficiencies 
found in large, strongly absorbing phytoplankton cells 
(Cannizzaro et al. 2008). 
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